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In the dense-neutrino region at 50-400 km above the neutrino sphere in a supernova, neutrino- 
neutrino interactions cause large flavor transformations. We study when the multi-angle nature of 
the neutrino trajectories leads to flavor decoherence between different angular modes. We consider 
a two-flavor mixing scenario between and another flavor v x and assume the usual hierarchy 
F Ve > F De > F„ x = Fv x for the number fluxes. We define e = (F Ve — Fp c )/(Fa c — Fv^.) as a measure 
for the deleptonization flux which is the one crucial parameter. The transition between the quasi 
single-angle behavior and multi-angle decoherence is abrupt as a function of e. For typical choices of 
other parameters, multi-angle decoherence is suppressed for e > 0.3, but a much smaller asymmetry 
suffices if the neutrino mass hierarchy is normal and the mixing angle small. The critical e depends 
logarithmically on the neutrino luminosity. In a realistic supernova scenario, the deleptonization 
flux is probably enough to suppress multi-angle decoherence. 



PACS numbers: 14.60.Pq, 97.60.Bw 

I. INTRODUCTION 

In the dense neutrino flux emerging from a supernova 
(SN) core, neutrino-neutrino refraction causes nonlinear 
flavor oscillation phenomena that are unlike anything 
produced by ordinary matter [l], 0, H, 0, H, H, 0, H, M, 
Hoi [nj. The crucial phenomenon is a collective mode 
of pair transformations of the form v e v e — > v x v x where 
x represents some suitable superposition of and v T . 
This pair-wise form of flavor transformation leaves the 
net flavor-lepton number flux unchanged. Even an ex- 
tremely small mixing angle is enough to trigger this ef- 
fect that is insensitive to the presence of ordinary matter 
unless there is a Mikheyev-Smirnov-Wolfenstein (MSW) 
resonance in the dense-neutrino region. 

Collective pair transformations require a large neutrino 
density and a pair excess of a given flavor. In typi- 
cal SN models one finds a hierarchy of number fluxes 
Fu e > F&e > Fi>x — Fv x - The first part of the hierar- 
chy is caused by the deleptonization of the collapsed core 
whereas the second is caused by the absence of charged- 
current interactions for neutrino species other than v e 
and v e . The neutrino fluxes streaming from a collapsed 
star thus provide a natural environment for a flavor pair 
excess. On the other hand, the SN core itself is character- 
ized by a large v e chemical potential that enhances the v e 
density and suppresses that of v e so that here pair trans- 
formations cannot occur. Likewise, the deleptonization 
burst immediately after core bounce has an excess of v e 
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and a depletion of v e [l2j , suggesting that it is unaffected 
by collective pair transformations. 

We illustrate collective pair conversions with a sim- 
ple example in Fig. [U assuming a typical SN neutrino 
luminosity to be quantified later. We show the evolu- 
tion of the z-components of the global flavor polarization 
vector P for neutrinos and P for antineutrinos, where 
initially P = |P| = 1 + e with e = 0.25 and P = |P| = 1. 
We have assumed a monochromatic spectrum, that all 
neutrinos are emitted at 45° relative to the radial direc- 
tion, the atmospheric Am 2 , a small vacuum mixing angle 
sin 26 = 10~ 3 to mimic the effect of ordinary matter, and 
an inverted mass hierarchy. For the normal hierarchy, no 
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FIG. 1: Schematic evolution of the z-components of the to- 
tal polarization vectors for neutrinos and antineutrinos in a 
SN caused by neutrino-neutrino interactions for the inverted- 
hierarchy example described in the text. 
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visible evolution takes place. 

Flavor oscillations do not change those parts of the fla- 
vor fluxes that are already equal, only the transformation 
of the excess v e flux over the v x flux is observable, and 
likewise for neutrinos. The polarization vectors only rep- 
resent this excess. Therefore, without loss of generality 
we may set F Vm = _Fp x = in our examples or equiva- 
lently, we may picture F Vc and F Pe to represent F„ e — Fp m 
and likewise for neutrinos. Our chosen parameters mean 
that at the neutrino sphere (R = 10 km) the excess of 
the v e flux over the v x flux is 25% larger than the ex- 
cess of the v e flux over the v x flux (e = 0.25). P z = +1 
then represents a pure u e excess flux, P z = represents 
equal excess fluxes of both flavors, and P z = — 1 a pure 
v x excess flux, and analogous for v e with f — > 1 + e. 

The main features of Fig. [T] are nicely explained after 
recognizing that the equations of motion can be brought 
into a form where they are equivalent to a gyroscopic pen- 
dulum 0, H| . The initial "plateau phase" corresponds to 
synchronized oscillations or, in the pendulum language, 
to a fast precession. We call the radius where this phase 
ends the synchronization radius r sync h. The decline with 
"wiggles" represents a nutation mode. The overall de- 
cline is caused by the dilution of the neutrino flux and 
their increasing collinearity with distance, corresponding 
to a decline of their effective interaction energy. 

One salient feature of Fig. [T] is that the v e flux com- 
pletely converts to v x , whereas the v e flux converts to 
v x only to the extent allowed by the conservation of 
P z — P z — e. This conservation is exact in the mass 
basis that approximately coincides with the interaction 
basis if the mixing angle is small. In other words, only 
v e v e pairs convert to v x v x pairs, whereas the unpaired v e 
excess remains in its original flavor Q. 

The current-current nature of the weak interaction 
causes the interaction energy to depend on (1 — cos 6*) for 
two trajectories with relative angle 0. Therefore, neutri- 
nos emitted in different directions from a SN core expe- 
rience different refractive effects % 6]. As a result, one 
would expect that their flavor content evolves differently, 
leading to kincmatical decoherence between different an- 
gular modes Q. Two of us have recently shown that 
this multi-angle decoherence is indeed unavoidable in a 
"symmetric gas" of equal densities of neutrinos and an- 
tineutrinos Q. Moreover, this effect is self-accelerating 
in that an infinitesimal anisotropy is enough to trigger an 
exponential run-away towards flavor equipartition, both 
for the normal and inverted hierarchies. 

In the SN context, however, it has been numerically 
observed that the evolution is much more similar to the 
single-angle (or the isotropic) case 0, 0] . The flux emit- 
ted by a SN is extremely anisotropic. If one assumes vv 
symmetry, flavor decoherence is swift and unavoidable. 
Therefore, the observed suppression of multi-angle deco- 
herence must be related to the v e v £ asymmetry that is 
generated by SN core deleptonization. 

To illustrate this point we show in Fig. [3 a few exam- 
ples along the lines of Fig. [H but now for multi-angle 



emission from the neutrino sphere that is again taken at 
10 km. We consider different values of an asymmetry 
parameter that we define as 

_ F{y e ) - F(v x ) _ F{v e ) - F(D e ) 

F(D e ) - F(D X ) F{v e ) - F(D X ) ' 1 1 

where we have used F(v x ) = F(y x ). As mentioned ear- 
lier, we can assume F(y x ) = F(D X ) = at the neutrino 
sphere without loss of generality. The left panels are for 
the normal hierarchy, the right panels for the inverted 
hierarchy. 

P z (r) — P z (r) — e is constant so that it is sufficient 
to show P z (r) alone. However, the length P — |P| is 
no longer preserved: Complete kincmatical decoherence 
among the angular modes would cause P = 0. On the 
other hand, if P = 1 remains fixed, this signifies that 
all modes evolve coherently with each other. We use P 
rather than P because the former measures what happens 
to the v e v e pairs, whereas the latter also includes the 
conserved v e excess. 

In the top row we use e = (symmetric case). The 
flavor content decoheres quickly as expected. Both the 
length and the z-components of P and P shrink to zero 
within about 20 meters of the nominal neutrino sphere. 

On the other extreme, we show in the bottom row the 
same for e = 0.25. In the normal hierarchy, nothing vis- 
ible happens, in analogy to the single-angle case. In the 
inverted hierarchy, the transformation is similar, but not 
identical, to the single-angle case. The nutations wash 
out quickly. Shortly after exiting from the synchroniza- 
tion phase, the length P shrinks a bit, but stays almost 
constant thereafter. Clearly, some sort of multi- angle ef- 
fect has happened as we will discuss further in Sec. IIIIl 
but multi-angle decoherence has certainly not occurred. 

In the two middle rows we show intermediate cases 
with e = 0.06 and 0.12, respectively. For the inverted hi- 
erarchy, these examples are qualitatively equivalent. The 
evolution is at first similar to the single-angle case and 
analogous to e = 0.25. The nutations are washed out and 
the length P shrinks a little bit after the synchroniza- 
tion radius. At some larger radius, however, something 
new happens in that P suddenly shrinks significantly, al- 
though not to zero, and there is a distinct feature in the 
evolution of the z-component. Now we obtain partial de- 
coherence. The final flavor content is very different from 
the single-angle case. 

In the normal hierarchy, and for e = 0.06, we obtain 
large decoherence that begins abruptly at some radius far 
beyond r sync h. For the larger asymmetry e = 0.12, the 
length P also shrinks, but closely tracks P z . As we will 
see, this case is somewhat like Phase II of the inverted- 
hierarchy case, i.e., a certain amount of shrinking of the 
length of P and thus a clear multi-angle effect, but no 
real decoherence. 

Depending on the deleptonization flux, here repre- 
sented by the asymmetry parameter e, the system be- 
haves very differently. In particular, for the inverted hi- 
erarchy it is striking that there are either two or three 
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FIG. 2: Radial evolution of P z in a schematic SN model as in Fig. [TJ but now for multi-angle neutrino emission at the neutrino 
sphere (R = 10 km). In addition we show the length P — |P| as a measure of kinematical coherence. Left: normal hierarchy. 
Right: inverted hierarchy. From top to bottom: e = 0, 0.06, 0.12 and 0.25, where e is defined in Eq. ([!]). 
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distinct phases. We always have the initial synchro- 
nized phase at large neutrino densities. Next, there is 
always the quasi single-angle pair-transformation phase 
at distances larger than r sync h. Just beyond this radius, 
the global polarization vectors quickly shrink by a small 
amount, but then stabilize immediately. Finally, if e is 
below some critical value, there is a sharp transition to a 
third phase where the different angular modes decohere 
significantly, but not completely. The practical outcome 
for the flavor fluxes emerging from the dense-neutrino 
region is very different depending on e. The transition 
between these regimes is abrupt, a small change of e is 
enough to cause one or the other form of behavior. 

While these phenomena call for an analytic quantita- 
tive understanding, we are here less ambitious, but more 
practical. We study numerically for which range of pa- 
rameters the different forms of behavior occur. Towards 
this goal we first set up, in Sec. UH our conventions, the 
equations of motion for a spherically symmetric system, 
and establish the connection between the parameters of 
our schematic model with those of a realistic SN scenario. 
In Sec. Hill we describe in more detail what happens in 
the different phases of evolution diagnosed in Fig. [2] and 
identify useful measures of decoherence. In Sec. lIVI we in- 
vestigate the role of our various model parameters in de- 
termining if the system kinematically decoheres. We dis- 
cuss our findings and conclude in Sec.fV] In Appendix El 
we derive the equations of motion adapted to spherical 
symmetry. 

II. SETUP OF THE PROBLEM 
A. Equations of motion 

To study the flavor evolution of the neutrino flux emit- 
ted by a SN core we solve numerically the equations of 
motion for the flavor-dependent number fluxes, assuming 
spherical symmetry We always work in a two-flavor sce- 
nario between is e and another flavor v x , characterized by 
the atmospheric Am 2 and by a vacuum mixing angle 9 
that is taken to represent the unknown I3-mixing angle. 

Our fundamental quantities are the flux matrices in 
flavor space J r that depend on the radial coordinate r 
(Appendix [A"|l . The diagonal entries represent the total 
neutrino number fluxes through a sphere of radius r. In 
the absence of oscillations, J r would not depend on the 
radius at all. The flux matrices are represented by po- 
larization vectors P r in the usual way, 

_ F(y e )+F(y x ) F{D e )~F{v x ) ^ 



- _ F(D e )+F(D x ) F{9 e )-F(y x ) ^ 

Jr — 2 ' 2 r ' ' 

where <r is the vector of Pauli matrices. Antineutrino 
quantities are always denoted with an overbar. The num- 
ber fluxes F(v) are understood at the neutrino sphere. In 
both equations the term proportional to the polarization 
vector is normalized to the antineutrino flux. As a conse- 
quence, at the neutrino sphere we have the normalization 

P=|P|=l + e and P = |P| = 1 . (3) 

In this way, we treat the excess flux from deleptonization 
as an adjustable parameter without affecting the baseline 
flux of antineutrinos. 

The diagonal part of the flux matrices is conserved and 
irrelevant for flavor oscillations. The polarization vector 
P,. only captures the difference between the flavor fluxes. 
For this reason we have defined the asymmetry e in terms 
of the flux differences. 

Multi-angle effects are at the focus of our study. We 
label different angular modes with 

u = sin 2 d R , (4) 

where is the zenith angle at the neutrino sphere r = R 
of a given mode relative to the radial direction. The pa- 
rameter u is fixed for every trajectory whereas the physi- 
cal zenith angle i? r at distance r varies. Therefore, using 
the local zenith angle to label the modes would compli- 
cate the equations. 

We will consider two generic angular distributions for 
the modes. In the multi-angle case we assume that the 
neutrino radiation field is "half isotropic" directly above 
the neutrino sphere, i.e., all outward moving modes are 
equally occupied as expected for blackbody emission. 
This implies (Appendix |A")) 

P«.r = dP r /du = const. (5) 

at r = R for < u < 1. Note that u = represents 
radial modes, u = I tangential ones. The other generic 
distribution is the single-angle case where all neutrinos 
are taken to be launched at 45° at the neutrino sphere 
so that u = 1/2 for all neutrinos. 

For a monochromatic energy distribution, the equa- 
tions of motion in spherical symmetry are (Appendix fA|) 



where the radial velocity of mode u at radius r is 

v u>r = y/l-uR 2 /r 2 . (7) 

Further, ui = \Am 2 /2E\ is the vacuum oscillation fre- 
quency, taken to be positive. B = (sin 26, 0, ± cos 29) 
where the mixing angle 9 is usually taken to be small. 
B z < corresponds to the normal hierarchy, B z > 
to the inverted hierarchy. L is a unit vector in the z- 
direction because we work in the interaction basis. The 
matter density is represented by 

A r = V2G F [n e - (r) - n e+ (r)] . (8) 

The strength of the neutrino-neutrino interaction is pa- 
rameterized by 

fi = V2G F (F«-F*) , (9) 

where the fluxes are taken at the neutrino sphere with 
radius R. 

The somewhat complicated structure of the equations 
arises from projecting the evolution of each mode on the 
radial direction. This is still very much simpler than 
following the evolution on every trajectory as a function 
of distance (or time) on that trajectory. We have here 
a closed set of differential equations that is not hard to 
solve numerically. 

We show in Appendix [X] that for r ^> R, the vacuum 
and matter oscillation terms take on the familiar plane- 
wave form because at large distances all neutrinos essen- 
tially move on radial trajectories. The neutrino-neutrino 
term falls off as r~ 4 , in agreement with the previous lit- 
erature. 



B. Schematic supernova model 

We always consider a two-flavor oscillation scenario 
driven by the atmospheric Am 2 — 1.9-3.0 x 10~ 3 eV 2 . 
Assuming (E v ) — 15 MeV, the oscillation frequency is 
oj = 0.3-0.5 km -1 . To be specific, we use 

/Am 2 \ , , , 

LJ = \-^E~/=°- 3]im (10) 

as a benchmark value in the monochromatic model. 

The total energy output of a SN is around 3 x 10 53 erg, 
corresponding to 0.5 x 10 53 erg in each of the six neutrino 
species if we assume approximate equipartition of the 



TABLE I: Default values for our model parameters. 



Parameter 


Standard value 


Definition 
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0.25 


Eq. © 


/' 


7 x 10 5 km" 1 


Eq. © 


u 


0.3 km" 1 


Eq. ([TO) 


sin 20 


10' 3 





emitted energy. If this energy is emitted over 10 s, the 
average luminosity per flavor would be 0.5 x 10 52 erg/s. 
However, at early times during the accretion phase, the 
luminosity in the v e flavor can exceed 3 x 10 53 erg/s [T^. 
As our baseline estimate we use 

fi = 7 X 10 5 km -1 (11) 

( L 9c \ 15 MeV / 10 km \ 2 

X {jE~) ~ JE~)J 10*2 erg/s \1T ) ' 

This is significantly larger than the assumptions of pre- 
vious studies 0,01 • Unless otherwise stated, we always 
use the benchmark values for the different parameters 
summarized in Table U 

In our calculations we always take the neutrino sphere 
at the radius R = 10 km. Of course, the physical 
neutrino sphere is not a well-defined concept. There- 
fore, the radius R simply represents the location where 
we fix the inner boundary condition. However, essen- 
tially nothing happens until the synchronization radius 
''synch 3> R because the in-medium mixing angle is ex- 
tremely small and both neutrinos and antineutrinos sim- 
ply precess around B. Therefore, as far as the vacuum 
and matter oscillation terms are concerned, it is almost 
irrelevant where we fix the inner boundary condition. 

Not so for the neutrino-neutrino term because we also 
fix the angular distribution at r = R. While the r~ 2 scal- 
ing from flux dilution is unaffected by the radius for the 
inner boundary condition, the "collinearity suppression" 
also scales as (R/r) 2 for r R. If we fix a half-isotropic 
distribution or a single angle of 45° at a larger radius R', 
the new inner boundary condition essentially amounts to 
fi — > fi' = n(R /R) 2 . In the early phase after bounce 
R' = 30 km could be more realistic, leading to a fi value 
almost an order of magnitude larger. Evidently, /i is a 
rather uncertain model parameter that can differ by or- 
ders of magnitude from our benchmark value. 

However, collective pair conversions only begin at 
''synch where fi is so small that synchronization ends. 
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Therefore, the main impact of a modified /i is to change 
r synch and thus to push the collective pair conversions to 
larger radii. The oscillations are synchronized if Q 



> 



(i - vr+ey 



(12) 



In our single-angle case we find from Eq. (|A20[) that the 
effective neutrino-neutrino interaction strength varies at 
large distances as 



R 

Therefore, the synchronization radius is 



(13) 



^synch 

R 



VT+e-l 



1/2 



(£) 



1/4 



2 W 



1/4 



(14) 



The second line assumes f«l. If we use our benchmark 
values uj = 0.3 km -1 , /i=7x 10 5 km -1 , R = 10 km 
and e = 0.25, we find r sync h = 95 km, corresponding 
well, for example, to Fig. [T] In any event, if fi is taken 
to be uncertain by two orders of magnitude, r sync h only 
changes by a factor of 3. 

The total electron lepton number emitted from a 
collapsed SN core is about 3 x 10 56 . On the other 
hand, assuming that each neutrino species carries away 
0.5 x 10 53 erg with an average energy of 15 MeV, the 
SN core emits about 2 x 10 57 neutrinos in each of the 
six species. In this simplified picture, the SN emits on 
average about 15% more v e than v e . However, in the 
oscillation context we need the excess of F Ve — F Vx rela- 
tive to the same quantity for antineutrinos as defined in 
Eq. (fTJ) . The true value of e thus depends sensitively on 
the detailed fluxes and spectra of the emitted neutrinos. 
The asymmetry parameter is large when the first hierar- 
chy in F Ve > F„ e > = F Vx is large and/or the second 
hierarchy is small. Even if F Px is as small as half of Fa e , 
the asymmetry e would be as large as 30%, even when 
F Ve exceeds F Pe by only 15%. 



C. Numerical multi-angle decoherence and the 
inner boundary condition 

One important and somewhat confusing complication 
of numerically solving the equations of motion is the phe- 
nomenon of numerical multi-angle decoherence. To inte- 
grate Eq. ^ one needs to work with a finite number of 
angular modes, equivalent to coarse-graining the phase 
space of the system. If the number of angular bins is 
chosen smaller than some critical number AT m i n , multi- 
angle decoherence occurs for r < r sync h, where physically 
it is not possible and does not occur for a fine-grained 
calculation. This phenomenon is shown, for example, in 
Fig. 3 of Ref. It is not caused by a lack of numeri- 
cal precision, but a result of the coarse-graining of phase 



space. A related phenomenon is recurrence as discussed 
in the context of multi-angle decoherence in Ref. [8|. 

In other words, a coarsely grained multi-angle system 
behaves differently than a finely grained one. A smaller 
mixing angle reduces AT m i n , a larger neutrino-neutrino 
interaction strength increases it. It should be possible to 
estimate N m - ln from first principles, but for the moment 
we need to rely on trial and error. 

Starting the integration at r — R is doubly punishing 
because the fast oscillations of individual modes caused 
by a large [i requires many radial steps for the numerical 
integration and avoiding numerical decoherence requires 
a large number of angular modes. On the other hand, 
in this region nothing but fast synchronized oscillations 
take place that have no physical effect if the mixing angle 
is small. Using a larger radius as a starting point for the 
integration avoids both problems and does not modify 
the overall flavor evolution at larger distances. 

From the physical perspective, the "neutrino sphere" 
is not a well-defined concept because different energy 
modes and different species decouple at different radii, 
and in any case, each individual neutrino scatters last at 
a different radius. If the exact inner boundary condition 
would matter, we would need to solve the full kinetic 
equations, including neutral-current and charged-current 
collisions. It is the beauty of the neutrino-neutrino fla- 
vor transformation problem that the real action begins 
at r sync h, significantly outside the neutrino sphere. Our 
approach of reducing the equations of motion to the re- 
fractive terms is only self consistent because the exact 
location of the inner boundary condition is irrelevant. 

In summary, the nominal neutrino sphere at R = 
10 km is nothing but a point of reference where we nor- 
malize the fluxes and fix the angular distribution. As 
a starting point for integration we typically use tq = 
0.75r sync h. A few hundred angular modes are then usu- 
ally enough to avoid numerical decoherence. 

We note, however, that the normal-hierarchy cases are 
more sensitive to both the number of angular bins and the 
starting radius for the integration. It can happen that a 
case that looks like the e = 0.12 example in Fig. [2J which 
shows a mild shrinking of the polarization vector, can 
become "more coherent" by choosing a smaller starting 
radius which then may also require a larger number of 
modes. For the normal hierarchy, the different multi- 
angle cases are less cleanly separated from each other 
than in the inverted hierarchy in that the transition is 
less abrupt as a function of e. 

When physical multi-angle decoherence occurs 
(e.g. the middle rows of Fig. [2]), a much larger number 
of modes is needed to provide reproducible results. 
However, we are here not interested in the exact final 
outcome, we are mostly interested in the range of 
parameters that lead to decoherence. Therefore, massive 
computation power is not needed for our study. 

For those cases where we include a non-trivial spec- 
trum of energies we also need energy bins. A distribution 
of energies does not lead to kinematical decoherence in 
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the context of collective neutrino oscillations [3] so that 
the number of energy bins is not a crucial parameter. 
Of course, to resolve the energy-de pen dent behavior and 
especially the spectral splits H, HllICj, Hp], a sufficiently 
fine-grained binning is required. It provides better reso- 
lution, but not a qualitatively different form of behavior. 



III. COHERENT EVOLUTION VS. 
DECOHERENCE 

A. Different forms of evolution 

Before investigating the conditions for decoherence 
among angular neutrino modes we first take a closer 
look at what happens in the different cases shown in 
Fig. O Considering first the quasi single-angle case with 
the asymmetry e = 0.25, some insight is gained by look- 
ing at the final state of the evolution at some large radius 
where the neutrino-neutrino effects have completely died 
out and all modes simply perform vacuum oscillations. 
In the left-hand panels of Fig. [3] we show the end state of 
500 polarization vectors, representing modes uniformly 
spaced in the angular coordinate u. In the upper panel 
we show the final state in the x-z-plane ("side view"), in 
the lower panel in the x-y-plane ( "top view" ) . 

Initially, all polarization vectors are aligned in the fla- 
vor direction. At the beginning of the pair transforma- 
tion phase at r sync h, some are peeled off, forming a spiral 
structure that is easily gleaned from the left panels of 
Fig. This structure continues to evolve almost as in 
the single-angle case, i.e., once established it moves al- 
most like a rigid body and eventually orients itself in the 
negative B-direction. Of course, it continues to rotate 
around the B direction even at large radii because of 
vacuum oscillations. 

The spiral structure is different depending on the mix- 
ing angle. We illustrate this in Fig. 2] where we show 
the top view in analogy to the lower-left panel of Fig. [3J 
for different choices of mixing angle. For a large sin 29, 
the polarization vectors stay close to each other. For a 
smaller sin 29, the spiral spreads over a larger solid angle 
and has more windings. We recall that a smaller sin 29 
also has the effect of causing a larger nutation depth of 
the flavor pendulum 

Now turn to the quasi decoherent case with e = 0.12. 
Initially the same happens, but at the "decoherence 
radius" the spiral structure dissolves almost instanta- 
neously. The polarization vectors enter a complicated 
structure as illustrated by the end state (central panels 
of Fig. [3]). Moreover, they are spread out all over the unit 
sphere, having both positive and negative z-components. 
This structure looks different for different choices of sin 29 
and e. However, once a sufficient number of polarization 
vectors is used, it is reproducible. For e = 0.06 the pic- 
ture would be qualitatively similar. 

Finally we show the fully symmetric case (e = 0) in 
the right-hand panels. Here decoherence is fast and com- 



plete. For a small mixing angle, all polarization vectors 
are confined to the x-z-planc. They distribute themselves 
on a circle in that plane. 

For the normal hierarchy, we show in Fig. [5] as an ex- 
plicit example the e = 0.12 case of Fig. [5] that showed 
a clear multi-angle effect without strong decoherence. 
Once more we find a spiral structure. Most polarization 
vectors remain oriented roughly in their original direc- 
tion, but in this case also with a tail of a few polarization 
vectors reversed. The quasi decoherent case (e = 0.06) 
and the symmetric system produce similar final pictures 
as the corresponding cases of the inverted hierarchy. 



B. Measures of decoherence 



Even in the quasi-decoherent cases the unit sphere is 
not uniformly filled with polarization vectors. Rather, 
in the mono-energetic case considered here, the occupied 
phase space is a one-dimensional subspace of the unit 
sphere. It is parameterized by the angular variable u 
and shows a clear line-like structure. This picture sug- 
gests to use the length of this line on the unit sphere 
as another global measure besides the length P to dis- 
criminate between different modes of evolution [8]. In a 
numerical run with discrete angular bins, this quantity 
is simply the sum of the angles between neighboring po- 
larization vectors. In Fig. [5] we show this quantity for 
the indicated values of e as a function of radius for our 
inverted-hierarchy examples. 

At the radius r sync h where the spiral forms, the length 
on the unit sphere quickly increases from to a value that 
is almost independent of e, but depends on the mixing 
angle. For smaller sin 29 it is larger, corresponding to the 
spiral having more windings as indicated earlier. Later, 
this length stays practically constant, reflecting that the 
spiral structure, once established, does not change much 
except tilting toward the negative B-direction and pro- 
cessing around it. 

When e is smaller than a critical value, at the "de- 
coherence radius" a sudden second growth phase shoots 
up from the plateau of these curves. For smaller e, the 
final length is longer, representing a more "phase-space 
filling" line on the unit sphere. 

Note, however, that for e close to zero, the line does 
not fill the unit sphere, but essentially stays in a narrow 
band. In the perfectly symmetric case, the motion of 
all polarization vectors is essentially confined to the x-z- 
plane, i.e., the polarization vectors distribute themselves 
over a great circle on the sphere as shown in the right 
panels of Fig. [3] 
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FIG. 3: Final location on the unit sphere of 500 antineutrino polarization vectors for our standard parameters and the inverted 
hierarchy. The top row is the "side view" (x-z-components), the bottom row the "top view" (rr-y-components). Left: quasi 
single-angle case (e = 0.25). Middle: decoherent case (e = 0.12). Right: symmetric system (e = 0). 
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FIG. 4: Same as Fig. [3] now only top views for quasi single-angle cases with the mixing angles sin 26 = 0.1, 10 -3 and 10 
from left to right. The middle panel is identical with the bottom left panel of Fig. [3] 



IV. ROLE OF MODEL PARAMETERS 
A. Ordinary matter 

We now explore how various model parameters influ- 
ence the behavior of the system. In the examples so far 
we have ignored matter because its effect is mainly to 
suppress the vacuum mixing angle. We here make this 
argument more precise. In Fig. [7] we show typical matter 
density profiles, expressed in terms of the matter oscil- 



lation frequency A(r), from numerical simulations of the 
Garching group for different times after collapse [l4[ • For 
comparison we also show /i(r) with fJ,(R) = 7 x 10 5 km -1 
and a radial variation in analogy to Eq. (|A20D . 

We observe that for the shown density profiles, the 
line uj intersects A(r) at a radius far exceeding the dense- 
neutrino region that lies within the radius where the fj,(r) 
profile intersects to. In other words, the H-resonance is far 
outside the region of interest except perhaps for very late 
times. Then, of course, the neutrino luminosity will be 
much smaller, i.e., the /i(r) curve would also shift down- 
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FIG. 5: Same as Fig. [3] here for the normal hierarchy and 
e = 0.12. 



ward and the dense-neutrino region would be limited to 
smaller radii. 

The true density profiles may be much lower, espe- 
cially at late times. This is even required for success- 
ful r-process nucleosynthesis. In this scenario an MSW 
resonance may take place within the dense-neutrino re- 
gion, a case that was the focus of previous numerical 
studies @|. However, we will always assume that the 



H-resonance is at larger radii and that neutrino-neutrino 
refraction and ordinary matter effects do not interfere. 

What is the impact of a large matter density in the 
region where neutrino-neutrino effects are important? In 
the previous literature it was recognized that a constant 
matter profile essentially reduces the effective mixing an- 
gle so that matter should have the same influence as a 
small vacuum mixing angle [J, [7| . We illustrate this point 
in Fig. [U with the evolution of P z for our usual case, but 
assuming now a large vacuum mixing angle sin 29 = 0.1. 
In the synchronization region one can now see oscilla- 
tions. We overlay this curve with P z , using the matter 
profile of Fig.[7]at t = 2 s. As expected, matter has the ef- 
fect of slightly delaying the onset of pair transformations 
and of increasing the depth of the nutation amplitude. 

Actually, in the inverted hierarchy, the value of sin 29 
is only crucial at the onset of the bipolar oscillations. 
Once the overall polarization vector is tilted away from 
B, the initial "misalignment" with B no longer matters. 
Therefore, what is crucial for the role of matter is only 
its density around the region where synchronization ends. 
The in-medium mixing angle at r sync h for the case shown 
in the top panel of Fig. [8] is sin20 mat tcr = 3.35 x 10~ 4 , 
assuming sm20 vac = 0.1. Using this value of # ma tter as 
a vacuum mixing angle instead of matter yields the re- 
sult shown in the bottom panel, again overlaid with the 
original vacuum case of sin 29 = 0.1. 

We conclude that indeed we can ignore matter entirely 
if we account for it schematically by a small vacuum mix- 
ing angle, at least in the inverted hierarchy. Moreover, 
the onset of collective pair transformations is only mildly 
changed by the choice of mixing angle. Its main impact 
is that it controls the depth of the nutation pattern. The 
exact matter profile is only important if it is so shallow 
that it causes an MSW resonance in the dense-neutrino 
region, a case that we do not investigate. 
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FIG. 6: Evolution of the length of the one-dimensional sub- 
space occupied by the polarization vectors for our standard in- 
verted hierarchy case, taking a series of different asymmetries 
e. The length grows to larger values for smaller asymmetries. 




FIG. 7: Typical matter density profiles from numerical sim- 
ulations of the Garching group at the indicated times after 
core bounce [lj. For comparison we also show our bench- 
mark value uj = 0.3 km -1 and /j.(r) for a typical case. 
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B. Mixing angle 

This discussion suggests that, at least for the in- 
verted hierarchy, the actual vacuum mixing angle does 
not strongly influence the issue of multi-angle decoher- 
ence because this effect happens when the global polar- 
ization vector is tilted far away from the B direction. On 
the other hand, we have already noted that the quasi- 
coherent spiral structure that forms just beyond the syn- 
chronization radius has more windings for a smaller mix- 
ing angle so that the system is not identical. 

To clarify the role of the mixing angle we have used our 
standard inverted-hierarchy case and have calculated the 
limiting asymmetry e for decoherence for a broad range 
of mixing angles. We show the limiting contours in the 
plane of e and sin 29 in Fig. [5] for both hierarchies, above 
which multiangle decoherence does not appear. 

We emphasize that the limiting e shown in Fig. [5] has 
a different meaning for the two hierarchies. As discussed 
earlier, in the inverted hierarchy, P shortens somewhat 
even in the quasi single-angle regime. Therefore, as a for- 
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FIG. 9: Limiting e for decoherence as a function of mixing 
angle for our standard example and both hierarchies. 



mal criterion for distinguishing the regions of coherence 
and decoherence we use that the final P has shortened 
to less than 0.85. The exact choice is irrelevant because 
the transition between the quasi-coherent and decoherent 
regimes is steep as a function of e. 

Conversely, in the normal hierarchy, P need not visibly 
shorten at all as illustrated by the example in the lower 
left panel of Fig. [21 Therefore, we here demand that P 
does not visibly shorten in such a picture. We construct 
the demarcation line by decreasing e in steps of 0.01 until 
the polarization vector for the first time shortens visibly. 
Finding this point requires a significant amount of man- 
ual iterations with a modified inner radius and number 
of angular bins to make sure the result does not depend 
on these numerical parameters. The error bars represent 
our confidence range for the true critical value. 

We conclude that for the inverted hierarchy, multi- 
angle decoherence is virtually independent from the value 
of sin2f5, except that for very large 9 a slightly smaller 
asymmetry is enough to suppress decoherence. Assuming 
the presence of ordinary matter, such large mixing angles 
seem irrelevant, except perhaps at late times. Either way, 
it is conservative to assume a small mixing angle and we 
will use sin 29 = 10 -3 as a default value. 

For the normal hierarchy we find a strong dependence 
of the critical e on log 10 (sin 29). For a smaller mixing 
angle it is easier to suppress decoherence. The normal hi- 
erarchy is very different from the inverted one in that for 
a small mixing angle, all polarization vectors stay closely 
aligned with the z-direction unless multi-angle decoher- 
ence takes place. Therefore, it is plausible that for a 
smaller mixing angle, decoherence effects arc delayed. 



FIG. 8: Evolution for our standard inverted-hierarchy case for 
a large vacuum mixing angle of sin 26 — 0.1 (blue/dotted lines 
in both panels), compared with different ways of suppressing 
the mixing angle (red/solid). Top: ordinary matter effect 
according to the profile at t — 2 s in Fig. [7] Bottom: small 
vacuum mixing angle of sin 2(9 = 3.35 x 10~ 4 . 



C. Energy distribution 

The neutrinos emitted from a SN core naturally have 
a broad energy distribution. In Ref. [8| it was noted that 
the energy distribution of neutrinos and antineutrinos is 
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FIG. 10: Final state at a large radius of the polarization vectors for our standard parameters in analogy to Fig. [3] The 
antineutrinos (red/light gray) are on the unit sphere, whereas the neutrinos (blue/dark gray) live on a sphere of radius 
1 + e = 1.25. Left: monochromatic multi-angle, the antineutrinos being identical with the left column of Fig. [3] Middle: 
Box-like energy spectrum and single angle. Right: Box-like energy spectrum and multi angle. In the lower right panel we do 
not show the antineutrinos. 



largely irrelevant for the question of decoherence as long 
as the oscillations exhibit self-maintained coherence [l5| . 
The multi-angle transition to decoherence typically oc- 
curs within the dense-neutrino region where the synchro- 
nization of energy modes remains strong. Therefore, we 
expect that multi-angle decoherence is not significantly 
affected by the neutrino spectrum. 

In order to compare a monochromatic system with one 
that has a broad energy distribution, the crucial quantity 
to keep fixed is not the average energy, but the average os- 
cillation frequency (u>) = (Am 2 /2E). If we assume that 
neutrinos and antineutrinos have equal distributions, it is 
straightforward to adjust, for example, the temperature 
of a thermal distribution such that (u) is identical to our 
monochromatic standard case luq = 0.3 km - . 

If we assume different distributions for neutrinos and 
antineutrinos, the equivalent loq is somewhat more subtle. 
Consider first two different monochromatic spectra for 
neutrinos with a fixed frequency uji , and one for antineu- 
trinos with a different frequency 0J2 ( "bichromatic sys- 
tem"). Following Ref. one can return to a monochro- 
matic situation by going into a reference frame that ro- 
tates around B with such a frequency that in vacuum 
P and P precess around B with equal frequencies wq, 



but in opposite directions. The rotation frequency for 
the corotating frame is oj c = (u>i — u>2)/2. Therefore, our 
bichromatic system behaves equivalently to a monochro- 
matic one with ojq = {lo\ + l>2)/2. It is trivial to show 
in numerical examples that the bichromatic system is in- 
deed equivalent to a monochromatic one with luq taken 
as the simple average of uj\ and 102- 

If we have different distributions for neutrinos and an- 
tineutrinos, we define the initial average frequencies by 



(15) 



and analogous for (ujp). The equivalent monochromatic 
frequency is then u> = i ((u> v ) + (u>p))- The initial dis- 
tribution P* can involve negative values if some part of 
the spectrum initially consists of v x and not v e . Such 
spectral cross-overs occur, for example, if one assumes 
thermal fluxes with equal luminosities but different tem- 
peratures. 

We have studied several numerical examples of quasi 
single-angle behavior and of multi-angle decoherence, 
taking different neutrino and antineutrino energy spec- 
tra, such as flat or thermal and with equal or different 
temperatures. We always found that the evolution of 
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the global polarization vectors is almost identical to the 
equivalent monochromatic cases. We never observed that 
a broad energy spectrum caused a significant deviation 
from the monochromatic behavior at those radii that are 
relevant for decoherence. 

Of course, a multi-energy system is qualitatively dif- 
ferent from a monochromatic one in that the final energy 
distribution shows a "spectral split" H H, H, Qil HH • ^ n a 
single-angle multi-energy system, this means the v e spec- 
trum is completely transformed to the v x flavor, whereas 
only the high-energy part of the v e spectrum is trans- 
formed, the low-energy part remaining in (or rather re- 
turning to) the original flavor. The energy E sp y lt of this 
sharp transition is fixed by lepton-number conservation 
in the sense that the neutrino-neutrino interactions only 
catalyzes the transformation of v e v e pairs. For various 
examples we find results in full agreement with the pre- 
vious literature % i, @, 0, EH • 

For sufficiently large asymmetries e where the multi- 
angle system evolves in the quasi single-angle mode, there 
is no significant modification of the spectral split so that 
it is not worthwhile to show any examples. In the de- 
coherent case, the final spectra naturally are very dif- 
ferent, but we have not explored such cases systemati- 
cally because multi-angle decoherence does not seem to 
be generic for realistic SN scenarios. 

To illustrate the modifications caused by an energy 
spectrum in a different way from the previous literature, 
we show in Fig. [10] the side and top views of the loca- 
tion of neutrino and antincutrino polarization vectors on 
the unit sphere in analogy to Fig. [3] for our standard pa- 
rameter values. In the left column we show the same 
monochromatic multi-angle case that we already showed 
in the left column of Fig. [3j with 500 modes. In addi- 
tion we include the neutrinos (blue/dark gray) that here 
live on a sphere of radius 1 + e = 1.25. The neutrinos 
form a spiral structure similar to the one of the antincu- 
trinos, but in the final state this structure cannot move 
to the negative B directions because of lepton number 
conservation. 

In the middle column we show a single-angle example 
with the same parameters, now using a box-like spectrum 
of oscillation frequencies where initially P* = (2ujq)~ 1 
and P* = (1 + e)(2w ) _1 for < w < 2w so that 
(w„) = (u}jj) = an d it is equivalent to the original 
monochromatic case. We now see that most of the an- 
tineutrinos have moved to the negative B direction as 
before, whereas the neutrinos populate both the posi- 
tive and negative B direction, representing the spectral 
split. The lack of full adiabaticity prevents the split 
from being complete, leaving some polarization vectors 
not fully aligned or anti-aligned with B. At large radii 
when the neutrino- neutrino interactions have died out, 
these modes precess with their different vacuum oscilla- 
tion frequencies so that they are found on a spiral locus 
extending from the "south pole" to the "north pole" that 
gets wound up further at larger radii. Note that here we 
have used 1000 energy modes in order to obtain a visi- 



ble population occupying these non-adiabatic final states. 
Still, only very few red dots (antineutrinos) are visible, 
the vast majority being at the south pole. Likewise for 
the neutrinos (blue dots), the spiral is populated only by 
a small fraction of the 1000 modes. In other words, the 
evolution is nearly adiabatic. 

Finally we combine a box-like energy spectrum and a 
multi-angle distribution (right panels). The antineutri- 
nos all cluster around the negative B direction and fill 
the "southern polar cap" more or less uniformly because 
at late times modes with different energies precess with 
different frequencies. The neutrinos populate both the 
northern and southern polar caps, representing the spec- 
tral split. At intermediate latitudes we find coherent spi- 
ral structures. They correspond to modes with different 
angles but equal uj so that even at late times they do not 
dissolve by differential precession. 



D. Effective interaction strength 

Besides the asymmetry e itself, the most uncertain 
model parameter is the effective neutrino-neutrino inter- 
action strength fj, as defined in Eq. ©. In Fig. [TT1 we 
show the demarcation lines between coherence and deco- 
herence for both hierarchies in the /i-e-plane, keeping all 
other parameters at their standard values. The contours 
are constructed as described in Sec. lIVBl The numerical 
contours are visually very well approximated by linear 
regressions of the form 



e m ~ 0.225 + 0.027 log 
e NH « 0.172 + 0.087 log 10 



'" l 10 6 km" 1 



10 6 km" 



(16) 



For the normal hierarchy, the linear regression would in- 
tersect e = within the range of investigated /^-values, 
but in reality turns over and saturates around e = 0.06. 



E. Vacuum oscillation frequency 

The average vacuum oscillation frequency lu depends 
on the atmospheric Am 2 that is quite well constrained, 
and a certain average of the neutrino energies. Our stan- 
dard value is to — 0.3 km -1 . If we increase this to 1 km -1 , 
the e-//-contour in Fig.[Tl]is essentially parallel-shifted to 
larger e by about 0.035 (inverted hierarchy). This range 
of to probably brackets the plausible possibilities so that 
the uncertainty of to does not strongly influence the prac- 
tical demarcation between the regimes. 

The normal hierarchy is more sensitive to to. In Fig.[T2l 
we show a contour for the coherence regime in the e-u> 
plane, assuming otherwise our standard parameter val- 
ues. Changing to from 0.3 to 1 km -1 increases the critical 
e by almost 0.15. 
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FIG. 11: Limiting e for decoherence as a function of the effec- 
tive neutrino-neutrino interaction strength for our standard 
parameters. The linear regressions are "visual fits" repre- 
sented by Eq. fig)). 
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FIG. 12: Limiting e for decoherence as a function of the vac- 
uum oscillation frequency uj for our standard parameters and 
the normal hierarchy. 



V. CONCLUSIONS 

The nonlinear neutrino transformations that occur in 
the dense-neutrino region of a SN show numerous novel 
features. It was noted that multi-angle effects play an 
important role in that the neutrino-neutrino interaction 
depends on the relative angles of the various trajecto- 
ries HQ. At the same time it was numerically observed 
that for a typical example the behavior was unexpectedly 
quite similar to the single-angle case 5, 6]. On the other 
hand it was analytically shown that a gas of equal den- 
sities of neutrinos and anti-neutrinos has a pronounced 
angular instability and kinematical decoherence between 
different angular modes in flavor space is fast, represent- 
ing a stable fixed point of the system Q . 

We have here not attempted to develop further ana- 



lytical insights, but have taken a practical approach and 
explored numerically the range of parameters where dif- 
ferent forms of behavior dominate in a realistic SN sce- 
nario. 

To this end we have first clarified that "multi-angle ef- 
fects" mean one of two clearly separated forms of behav- 
ior. The flavor content of the system can evolve in a quasi 
single-angle form. On the level of the polarization vectors 
this means that they fill only a restricted volume of the 
available phase space and maintain a coherent structure. 
On the other hand, nearly complete flavor equilibrium 
can arise where the available phase space is more or less 
uniformly filled. 

For realistic assumptions about supernova and neu- 
trino parameters, the switch between these modes of evo- 
lution is set by the degree of asymmetry between the neu- 
trino and antineutrino fluxes. While this asymmetry is 
caused by the deleptonization flux, the crucial parameter 
e is the asymmetry between F(v e ) — F(v x ) and the corre- 
sponding antineutrino quantity as defined in Eq. JT]) be- 
cause for flavor oscillations the part of the density matrix 
that is proportional to F{v e ) + F(v x ) drops out. While in 
a realistic SN on average F(v e ) is about 15% larger than 
F(9 e ), the asymmetry parameter as defined in Eq. ([1]) is 
typically much larger. 

The critical value of e that is enough to suppress deco- 
herence depends on the type of neutrino mass hierarchy, 
the average energies, luminosities, and on the mixing an- 
gle. We have found that for e > 0.3, decoherence is sup- 
pressed for the entire range of plausible parameters, but 
a value smaller than 0.1 may be enough, depending on 
the combination of other parameters. 

We conclude that the quasi single-angle behavior may 
well be typical for realistic SN conditions, i.e., that the 
deleptonization flux is enough to suppress multi-angle de- 
coherence. To substantiate this conclusion one should 
analyze the output of numerical simulations in terms of 
our model parameters. Besides the flavor-dependent lu- 
minosities and average energies, one needs the angular 
distribution of the neutrino radiation field at some ra- 
dius where collisions are no longer important. 

If our conclusion holds up in the light of realistic SN 
simulations, a practical understanding of the effect of self- 
induced neutrino flavor transformations quickly comes 
into reach. In the normal mass hierarchy, nothing new 
would happen on a macroscopic scale. In the inverted 
hierarchy, the final effect would be a conversion of v e v e 
pairs and a split in the v e spectrum. These phenomena 
are only mildly affected by multi-angle effects as long as 
we are in the quasi single-angle regime. 

If at late times the matter density profile contracts 
enough that an MSW effect occurs in the dense-neutrino 
region, the situation becomes more complicated as the 
neutrino-neutrino and ordinary matter effects interfere 
and produce a richer structure of spectral modifica- 
tions 5, 6]. Even then, numerical simulations are much 
simpler if multi-angle decoherence is suppressed. 

It is not obvious how e evolves at late times. The delep- 
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tonization of the core is probably faster than the cooling 
so that one may think that e becomes smaller. On the 
other hand, the v e can essentially only interact via neu- 
tral current reactions and their flux and energy distribu- 
tion should, therefore, become very similar to the ones 
of v x and D x . Therefore, it is not obvious if at late times 
the initial flux difference F Vii — or F Pc — i 7 ^ decreases 
more quickly. We also note that there can be a cross- 
over in the sense that at late times the flux hierarchy can 
become F Ux — F^ > F Vc > F Pc as in Ref. meaning 
that we would have a pair excess flux of v x v x instead of 
a v e v e excess. 

We have refrained from an interpretation of our numer- 
ical findings because we do not have developed a theory 
of kinematical decoherence for a system that is asym- 
metric between neutrinos and antineutrinos and where 
the effective interaction strength varies as a function of 
time (or here of radius). The absence of multi-angle de- 
coherence seems to follow from a lack of time for it to 
develop. One can interpret our results such that a more 
adiabatic decrease of the neutrino-neutrino interaction 
strength requires a larger asymmetry to suppress deco- 
herence. The different length scales of the problem seem 
to conspire such that the evolution is adiabatic in that 
sharp spectral splits develop, but not so adiabatic that 
kinematical decoherence would be typical. An analytic 
understanding of this conspiracy remains to be found. 

Our results suggest that signatures of collective flavor 
transformations are not erased by multi-angle decoher- 
ence and will survive to the surface, modulated by the 
usual MS W flavor conversions [16| . The survival of ob- 
servable signatures then also depends on the density fluc- 
tuations of the ordinary medium that can be a source of 
kinematical flavor decoherence [13, Gj| ■ 

All authors in this field have relied on the simplifying 
assumption of either homogeneity or exact spherical sym- 
metry to make the equations numerically tractable. The 
neutrino emission from a real SN is influenced by den- 
sity and temperature fluctuations of the medium in the 
region where neutrinos decouple. Likewise, the neutrino 
fluxes emitted from the accretion tori of coalescing neu- 
tron stars, the likely engines of short gamma ray bursts, 
have fewer symmetries than assumed here. It remains to 
be investigated if systems with more general geometries 
behave qualitatively similar to the spherically symmetric 
case or if deviations from spherical symmetry can provide 
a new source of kinematical decoherence. 
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APPENDIX A: EQUATIONS OF MOTION 

1. Temporal evolution 

A homogeneous ensemble of unmixed neutrinos is rep- 
resented by the occupation numbers f p = (a p a p ) for each 
momentum mode p, where op and a p are the relevant cre- 
ation and annihilation operators and (. . .) is the expecta- 
tion value. A corresponding expression can be defined for 
the antineutrinos, / p — (a p a p ), where overbarred quanti- 
ties always refer to antiparticles. In a multiflavor system 
of mixed neutrinos, the occupation numbers are gener- 
alised to density matrices in flavor space 0, H^, HH 

{Q P )ij = (a\aj) p and [g p ) l3 = (at a,) p . (Al) 

The reversed order of the flavor indices i and j in the 
right-hand side for antineutrinos assures that g p and g p 
transform identically under a flavor transformation. 

Flavor oscillations of an ensemble of neutrinos and an- 
tineutrinos are described by (l9l . l20l . |2~H 

idtgp = [Hp, Q p ] and id t g p = [H p , g p ] , (A2) 

where [•, •] is a commutator. The "Hamiltonian" for each 
mode is 

Hp = fl p + XL + V2G F J (ft, - ft,) (1 - v q ■ v p ) , 

(A3) 

where Gp is the Fermi constant. The matrix of vacuum 
oscillation frequencies for relativistic neutrinos is in the 
mass basis f2 p = diag(mf , m|,mf)/2p with p = |p|. The 
matter effect is represented by A = ^/2Gp(n e - — n e +) 
and L = diag(l,0,0), given here in the weak interaction 
basis. We ignore the possible presence of other charged- 
lepton flavors. The Hamiltonian for antineutrinos H p is 
the same with VL p — > — fi p , i.e., in vacuum antineutrinos 
oscillate "the other way round." 

The factor (1 — v q ■ v p ) = (1 — cos# pq ) represents 
the current-current nature of the weak interaction where 
v p = p/p is the velocity. The angular term averages to 
zero if the gas is isotropic. We ignore a possible net flux 
of charged leptons lest the ordinary matter effect also 
involves an angular factor. 

If the system is axially symmetric relative to some di- 
rection, the angular factor simplifies after an azimuthal 
integration to [5J, [8| 

(1 - v q ■ v p ) -» (1 - U qUp ) , (A4) 

where the velocities are along the symmetry axis. 
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2. Spatial evolution in spherical symmetry 

Instead of a homogeneous system that evolves in 
time we consider a stationary system that evolves in 
space. The occupation numbers become Wigner func- 
tions, which depend both on spatial coordinates and on 
momenta, but there is no conceptual problem as long as 
we consider spatial variations that are slow on the scale 
of the inverse neutrino momenta. 

Since multi-angle effects are at the focus of our prob- 
lem, we cannot reduce the equations to plane waves mov- 
ing in the same direction. Motivated by the SN applica- 
tion, however, we can take advantage of global spherical 
symmetry, implying that the ensemble is represented by 
matrices that depend on a radial coordinate r, the zenith 
angle relative to the radial direction, and the energy E 
which in the relativistic limit is identical with p = |p|. 

We ignore gravitational deflection near the SN core and 
assume that neutrinos move on straight lines after being 
launched at a radius R that we call the neutrino sphere. 
Consider a neutrino that was launched at an angle flu 
relative to the radial direction. Its radial velocity is 



V R = COS $ } 



(A5) 



At r > R the trajectory's angle relative to the radial 
direction is implied by simple geometry to be Q (see e.g. 
their Fig. 1) 



R sin i?r ; = r sin i? r 
Therefore, the radial velocity at r is 



cos d r 



where we have introduced 



R? 
1 u 



u = 1 



ti n = sin 



0«. 



(A6) 



(A7) 



(A8) 



It is convenient to label the angular modes with u. The 
physical zenith angles change with distance so that the 
equations would be more complicated. 

The density matrices g PtU ,r & re not especially useful 
to describe a spherically symmetric system because they 
vary with r even in the absence of oscillations. (Note 
that we often write the dependence of a quantity on a 
variable as an subscript.) A quantity that is conserved 
in the absence of oscillations is the total flux matrix 



J r ~ R? 



d 3 p 

(2tt) 3 



Qp.r ^p,r • 



(A9) 



To express the integral in co-moving variables we observe 
that d 3 p in spherical coordinates is p 2 dpdipdcosi} r and 
that Eq. JA7| implies 



d cos d r 



du 



1 



2v„ 



R*_ 



(A10) 



Therefore, we finally define the differential flux matrices 



P Qp,-, 



2(2tt) 



(All) 



where we have used J dip — 2ir for axial symmetry. The 
normalization is 



du 



(A12) 



In the absence of oscillations the total and differential 



fluxes are conserved, 9 r J r — and <9 r J 



r^p,u,r 



0. 



To include oscillations, we note that the radial velocity 
along a neutrino trajectory is v u-r = dr u /dt = cosi? UjT .. 
Therefore, if we wish to express the temporal evolution of 
the neutrino density matrix along its trajectory in terms 
of an evolution expressed in terms of the radial coordi- 
nate r, we substitute d t — > v ur d r in Eq. (|A2|) so that 



[H 



(A13) 



and analogous for antineutrinos. In other words, we 
project the evolution along a given trajectory to an evo- 
lution along the radial direction. For vacuum oscillations 
this has the effect of "compressing" the oscillation pat- 
tern for non-radial modes, i.e., even for monochromatic 
neutrinos, the effective vacuum oscillation frequency de- 
pends on both r and u. 

The vacuum-oscillation and ordinary-matter contribu- 
tions to V\p, u ,r were given in Eq. (|A3|) . whereas the self- 
term must be made explicit. To this end we introduce 
the matrix of number densities 



(A14) 



and its integral as 



f 1 f°° f°° f°° Jp«r 

H r = du dp N PiU r = / du dp ' 

JO JO Jo Jo v u 

(A15) 

Collecting all terms and taking advantage of Eq. (|A4|) for 
axial symmetry, we find 



id r Sp, u ,r — + \flp, Np jtl . r ] + \ r [L, Np iU>r ] + \/2Gf ( [Nr — N r , Np Ujr ] — [j r — J r , Jp,u,r] ^ 5 

R 2 / \ 

i<9rJp,u,r = -[fip, Np,«,r] +A r [L, Np tUtr ] + V2G F (jNr- N r , N p ,u,r] - [Jr~ Jr, Jp,«,r]J : ( A16 ) 
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where the electron density's radial variation is included 
in A r . 



3. Angular emission characteristics 

In a numerical simulation we need to specify the fluxes 
at the neutrino sphere r = R. For our usual multi- 
angle simulations we assume that the neutrino radia- 
tion field is "half isotropic" directly above the neutrino 
sphere, i.e., that all outward- moving angular modes are 
equally occupied as behooves a thermal radiation field. 
Therefore, the occupation numbers are distributed as 
dn/dcos$i?, = const., implying that the radial fluxes are 
distributed as dj/dcosi?^ = I'ijdn/dcostffl' oc cosi9r be- 
cause vr = cos^r- Expressed in the angular variable 
u this implies dj/du = const, because of Eq. (|A10[) . In 
other words, a blackbody radiation field at the neutrino 
sphere implies that 



J u = const. 



(A17) 



in the interval < u < 1. 

To avoid multi-angle effects one may sometimes wish 
to use a single angular bin. To represent a uniform J. u 
distribution, the natural choice is u = 1/2, corresponding 
to a launch angle da = 45°. Our numerical single-angle 
examples always correspond to this choice in an otherwise 
unchanged numerical code. 

In this case the radial velocity of all neutrinos as a 
function of radius is 



1 



2r 2 



(A18) 



For a monochromatic spectrum, the remaining flavor ma- 
trices are simply the total J r (corresponding to the single 
u = 1/2) and N r = } r /v r . Ignoring the trivial ordinary 
matter term, the equations of motion are 



[«,Jr] 



R 2 



1 



1 [J, 



(A19) 



and analogous for the antineutrinos. The coefficient of 
the neutrino-neutrino term is explicitly 



1 



R 2 /i 



(A20) 



At the neutrino sphere it is equal to ^/2G Fl whereas at 
large distances it is (v2Gt/2) R 4 /r 4 . As observed in the 
previous literature, the neutrino-neutrino term dies out 
at large distances as r~ 4 . 

One can define a "single-angle case" somewhat differ- 
ently. Assuming all angular modes evolve coherently, we 
can integrate the equations of motion over du and 

study the evolution of the quantities J p . r = J Q duJ p ^ u ^ r . 
To write the equations in a compact form we introduce 
the notation 



-4= -J- f°d P f du ] -^ 

v * ■'r JQ JQ v u,r 



(A21) 



The full equation of motion Eq. (|A16[) for neutrinos be- 
comes 



R 2 ( 1 



r 2 \(v*y 



1 J [ Jr Jr i Jp,r] 



(A22) 



and analogous for antineutrinos with fl p — > —Q p . 

At large distances we have 1/v* — 1 + \{R/r) 2 (u) 
where (it) is the average of u at emission. For the vacuum 
and matter terms, we only need the leading terms so that 
we recover the familiar plane-wave form of the equations 
of motion. The coefficient of the neutrino- neutrino term, 
on the other hand, becomes 



V2G F ^ («) 
r 4 



(A23) 



Both for half-isotropic emission and for our single-angle 
case we have (u) = 5, in agreement with Eq. (|A20[) . 
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